Abstract. We have performed Monte Carlo simulations 
of Compton upscattering of low-energy photons in an ac- 
cretion disk around a Schwarzschild black hole. The pho- 
tons gain energy from the rotational motion of the elec- 
trons in the disk. The upscattering occurs near the black 
hole horizon, where the flow velocity of the electrons ap- 
proaches the speed of light. We show that this type of 
bulk-flow Comptonization can produce power-law X-ray 
spectra similar to the ones observed in black-hole X-ray 
transients in the high/soft state, i.e., a soft bump dom- 
. inating the spectrum below ~ 10 keV and a power-law 
C^) tail with photon index in the range 2-3. In order to re- 
P\l ■ produce the observed hard to soft flux ratio the disk has 
$_i ' to have vertical optical depth above ~ 3 at the last stable 
Qh, orbit. We conclude that the power-law component of the 
' high/soft state of black-hole transients may be due to an 
^-j- , intrinsically cool disk extending all the way to the hole, 
C*\J ' without a separate hot plasma component. 
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1. Introduction 

The X-ray continuum spectrum of black-hole binaries is 
characterized by an ultrasoft component and a power-law 
hard tail. The ultrasoft component is interpreted as com- 
ing from an optically thick accretion disk. The temper- 
ature of the disk is high (<~ 1 keV) near the black hole 
and decreases outwards. Thus, the spectrum emitted by 
the disk is a multi-temperature black-body spectrum. The 
power-law tail is thought to originate from Comptoniza- 
tion of soft photons scattering off very energetic electrons. 
The amplitude of this component has been seen to corre- 
late with the 1 — 10 keV flux. During the low (also called 
hard) state, the ultrasoft component is very weak or ab- 
sent, whereas in the high (soft) state it dominates the 
spectrum, especially at energies <10 keV. 

Two types of Comptonization processes have been put 
forward to explain the power-law component: thermal and 
bulk-flow Comptonization. Thermal Comptonization re- 
quires the presence of a hot and rarefied cloud of electrons 
or corona (see, e.g., the review by Poutanen 1998 and ref- 
erences therein), where the soft photons get upscattered 
due to collisions with Maxwellian distributed 50-100 keV 
electrons. In the bulk-flow Comptonization scenario (see, 
e.g., the review by Kylafis & Reig 1999 and Papathanas- 
siou & Psaltis 2001), photons gain energy from collisions 
with radially infalling electrons. Titarchuk et al. (1997) as- 
sessed the relative importance of these two processes and 
concluded that the bulk motion is more efficient in upscat- 
tering photons than thermal Comptonization in spherical 
accretion, provided that the electron temperature in the 
radial flow is less than a few keV. 

Although these Comptonization models have been suc- 
cesful in describing the spectra of black-hole candidates, 
there are still a number of unresolved issues. A radial in- 
flow of plasma has been considered for mathematical sim- 
plicity, but how this radial inflow is formed is unclear. 
Since the disk is likely to have a magnetic field generated 
by the shear flow, it may well have a corona heated by 
magnetic dissipation (Galeev et al. 1979). Comptonization 
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Fig. 1. Variation of the vertical optical depth with radius 
in the inner disk (normalized to the optical depth to in 
the outer parts of the disk), for the model used in the 
scattering calculations. 



in this hot corona is one of the currently popular models 
for the hard spectral components in X-ray binaries, the 
other being two-temperature flows or 'ADAFs'. There is 
little independent evidence for cither of these possibilities, 
however, apart from the interpretation of the X-ray spec- 
trum they offer. 

On the other hand, general agreement exists on the 
presence of the accretion disk itself around a black hole in 
a mass transfering binary system. In this work we investi- 
gate the Comptonization process in a cold accretion disk 
around a black hole by Monte Carlo simulations. Soft pho- 
tons produced inside or outside the disk get upscattered 
by collisions with fast-moving electrons, which follow cir- 
cular trajectories outside the last stable orbit and ballistic 
trajectories inside it. We refer to this process by the name 
'orbital Comptonization'. It has been considered before 
by Hanawa (1990) who applied it to the boundary layer 
between a disk and a neutron star surface. 
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1.1. Thermal vs. bulk- flow Comptonization 

In a geometrically thin disk, the scattering photons escape 
easily before they have a chance to pick up energy from 
the orbital motion. Bulk-flow Comptonization is therefore 
likely to work better in thicker disks. Such disks are also 
hotter however, so that thermal Comptonization will also 
take place. 

To compare the relative contributions of the two pro- 
cesses, note that the scattering photons drift from their 
point of origin by random walk. Before they escape 
through the disk surface, they therefore sample the or- 
bital velocity field over a radial distance which is of the 
same order as the disk thickness. The typical velocity 
difference encountered by a scattering photon is thus of 
the order v — £IH ~ c s , where c s is the sound speed in 
the disk and O is the angular velocity. In the absence of 
radiation pressure, c s ~ (feT/mp) 1 / 2 . In thermal equi- 
librium, the thermal velocity of the scattering electrons 
is Csfmp/iTic) 1 ' 2 ~ (m p /m e ) 1/ ' 2 v ^> v. It would thus 
seem that thermal Comptonization wins over bulk-flow 
Comptonization by the orbital motion. For cool Shakura- 
Sunyaev disks this would indeed be the case. However, 
the temperatures in such disks, of the order 1 keV, are 
not sufficient for significant Comptonization anyway. 

More interesting are the geometrically thick accretion 
flows that are possible in the inner regions of disks. Here 
the conditions are different from those in standard thin 
disks. In the optically thin ADAF type flows, thermal equi- 
librium between protons and electrons does not hold, and 
the scattering electrons are much colder than expected 
from the thickness of the disk. The electron temperatures 
are near 100 keV so that the electron velocities are around 
0.5c. This is of the same order as the range of orbital veloc- 
ities sampled by a photon escaping through a disk of thick- 
ness H/r ~ 1 near the last stable orbit. In the optically 
thin type of ADAF, orbital Comptonization can therefore 
be as important as thermal Comptonization, though it 
may not dominate. 

Another type of radiatively inefficient accretion flow 
into a black hole takes place when the optical depth of 
the flow is large, and the disk is puffed up to aspect ratios 
H/r ~ 1 by radiation pressure. The temperature in such 
radiatively supported disks is relatively low, so that ther- 
mal Comptonization is not very effective, but the escaping 
photons will have sampled a large range of orbital veloci- 
ties. The conditions in optically thick radiation supported 
flows such as may be present at high accretion rates onto 
black holes are thus especially relevant for orbital Comp- 
tonization. 

2. The model 

2.1. Densities and velocity field 

In an accretion disk around a non-rotating black hole of 
mass M matter follows nearly circular orbits only up to 



the last stable orbit at a distance of n as t = 6r g , where 
r g = GM/c 2 is the gravitational radius. At the black hole 
horizon rj, = 2r g , the speed of the particles is c. 

The details of bulk Comptonization depend on the 
velocity field of the accretion flow, as well as its optical 
depth and geometrical thickness as a function of the dis- 
tance from the hole. For the high accretion rates we have 
in mind, the flow has a finite thinkness of the order r g , 
depending on the details of disk viscosity and radiative 
transfer. Instead of a more detailed quantitative model, 
we approximate the velocity field by purely circular mo- 
tion outside the last stable orbit, and a ballistic spiral-in 
inside r = n ast , as follows: 

For r > n ast , matter in the disk is taken to be in 
circular motion. In Schwarzschild coordinates, the orbital 
frequency as seen by a distant inertial observer is then Q = 
(GM/r 3 ) 1 / 2 [Misner, Thorne, & Wheeler 1973, (hereafter 
MTW) p671], the same as in Newtonian gravity. At r = 
riast and in geometric units, the energy of the particles per 
unit mass is E — y/8/9 and their angular momentum per 
unit mass L divided by the black- hole mass M is L/M = 
VT2 (MTW p662). 

To find the velocity of the g QlS clS el function of r as it 
falls ballistically into the black hole, we use the expressions 
(Shapiro & Teukolsky 1983 pp341-342) 

mo - if^y^ij 

We assume that just inside Hast, E = y/8/9 and L/M = 
VT2 — e, where e <C yl2 is a measure of the radial drift 
speed in the accretion disk just outside the last stable 
orbit. The exact value of e has no effect on the calculations. 

To compute the density and the vertical optical thick- 
ness of the flow, we assume a constant thickness H = 2r g 
(as seen by a distant observer). From the continuity equa- 
tion and the assumed velocity field, the vertical optical 
depth can then be computed from the radial velocity 
(Eq. |). Outside ri as t, where we have assumed circular 
orbits, we take a constant density. The resulting vertical 
optical depth of the disk as a function of distance, T z (r), 
is shown in Fig. |l|. 

Bulk-flow Comptonization that produces hard X-rays 
is effective only near the black-hole horizon, where the 
speed of the electrons is high. Thus, only those photons 
that find themselves in the inner part of the disk, and 
manage to escape eventually, will have their energy sig- 
nificantly increased. It is these photons that produce the 
high-energy power-law spectra. 

3. Radiative transfer 

Photons are generated throughout the flow, diffuse over 
some distance and leave it at some other point. The radi- 
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Fig. 2. Emergent spectral energy distributions as func- 
tions of position r/r g of the source of soft photons 
(monochromatic at 1 keV, in the midplane of the disk), 
and t , a measure of the vertical optical thickness (see text 
and Figure 1). To = (dot-dashed), 0.1 (short-dashed), 1 
(long-dashed), 3 (dotted), and 10 (solid). The three pan- 
els are for r = 3r g (top), r = 6r g (middle) and r = 10r g 
(bottom). 

ation leaving the disk at any point of its surface is thus an 
integral over the contributions from all source positions. 
In order to highlight the importance of different source 
positions (in particular, of their distance to the hole), we 
show here results for fixed source positions. The net spec- 
trum can be obtained from this by integration, if needed. 
The emission of photons at the source is assumed to be 
isotropic in a comoving frame, as it would be for photons 
generated by bremsstrahlung, for example. 

As a photon travels in the disk, it experiences Comp- 
ton scatterings with the orbiting electrons. If the optical 
depth is small (<1), the majority of the input photons 
escape unscattered. Those that are directed towards the 
black hole get absorbed. Those that scatter a few times 
have on average their energy significantly increased. If the 
optical depth is moderate, then those photons that are 
not absorbed by the black hole random-walk through the 
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Fig. 3. Same as in Fig. ^, but with a blackbody source 
(dot-dashed line) of soft photons of temperature kT = 0.5 
keV instead of a monochromatic source. 



medium prior to escape, and gain energy from the bulk 
motion of the electrons. 

We make some simplifications in dealing with the gen- 
eral relativistic aspects of the problem. The effect of cur- 
vature on the path of the photon between scatterings is 
neglected, that is, we approximate it as a straight-line as 
seen by a distant observer. Since the distances traveled 
between scatterings are small, this is not a bad approxi- 
mation. For a quantification of the small error made with 
the use of straight-line photon trajectories the reader is 
referred to Papathanassiou & Psaltis (2001). The effect of 
curvature on the geometrical distances is also neglected 
since the actual physical distance between scatterings is 
not important for the scattering process, only the physi- 
cal conditions at the scattering centers. Finally, we leave 
out the self-illumination of the disk surface by other parts 
of the disk (through gravitational light bending, for exam- 
pie). 

We do take into account the important effect of the 
changing gravitational redshift along the photon trajec- 
tory. The energy E(ri) of a photon, as measured by a lo- 
cally inertial observer at radius r\ , is related to the energy 
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Fig. 4. Dependence of the spectrum on disk thickness, 
H = 0.3 r g (dot-dashed line), 1 r g (dashed line) and 5 r g 
(solid line). The input source of photons is monochromatic 
at 1 keV and is located at 6 r g . The optical depth was 
taken to be 3 
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Fig. 5. Photon-number spectral index in the range 20-100 
keV as a function of vertical optical depth tq = r z (riast) f° r 
the spectra shown in Fig. || top (circles), middle (squares) 
and bottom (triangles). 



of the same photon Efa) at radius r2, as measured by a 
locally inertial observer through the difference in gravita- 
tional redshift (MTW p659) 

E( ri )y/l- 2r g /n = E(r 2 )y/l- 2r s /r 2 (3) 
3.1. Monte Carlo implementation 

The basic principle of the Monte Carlo technique, as it is 
applied to the Comptonization process (Cashwell & Ev- 
erett 1959; Pozdnyakov, Sobol, & Sunyaev 1983), consists 
of following the life histories of a large number of photons 
from the moment each photon is emitted from a given 
source until it leaves the scattering medium. 

The photon is characterized by 4 parameters (position, 
direction, energy and weight) which are updated at each 
scattering, and determine the spectrum of the radiation 
emerging from the scattering medium. The procedure can 
be divided into the following steps: 

1. The disk is divided into concentric cylindrical zones. 
Each zone is characterized by a mean velocity and a 
mean density. The number of zones is chosen so as to 
have a fine zoning of electron velocities in the inner 
parts of the disk. 

2. We calculate the optical depth to electron scatter- 
ing T(r,n,E) from the scattering position r to the 
boundary of the disk along the direction of the photon 
n. Then, a portion We~ T of the photon's weight es- 
capes and it is recorded, while the remaining weight 
W = W(l — e~ T ) is scattered between r and the 
boundary. 

3. We determine the new scattering position, the energy 
of the photon after scattering and its new direction. 
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Fig. 6. Hard to soft flux ratio as a function of the vertical 
optical depth tq — r z (ri as t) and source position: r = 3r g 



(circles), 



6r g (squares). The hard and soft bands are 



20-200 keV and 1-5 keV, respectively. The solid lines en- 
close the region compatible with observations. 

The process continues with the next scattering until 
the photon's weight falls below a certain value (10 -8 in 
our runs), and it is repeated for N photons, where TV is a 
large number (typically 10 7 ). 

4. Results and discussion 

We consider a point source of soft photons in the midplane 
of the accretion disk. In the local rest frame of the flow, 
the source emits monochromatic photons of energy Eq = 1 
keV isotropically. In the frame of a local inertial observer 
the distribution of emitted photons is not isotropic but 
forward peaked (in the direction of the flow) and the initial 
energy of the photons is Doppler shifted and depends on 
their initial direction. 
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The resulting spectral energy distributions E 2 f(E), 
where f(E) is the photon-number spectrum, are shown 
in Figs. H - [| Figure || shows results for three distances of 
the soft photon source from the hole, r/r g — 3, 6 and 10, 
and for five values of the vertical optical depth at the last 
stable orbit. All five spectra in each panel are normalized 
so that the integral of f(E) from zero to infinity is equal 
to 1. 

The emergent spectra in Fig. || consist of a peak, whose 
width increases as the source of photons approaches the 
black hole, and a high-energy tail that is approximately 
a power law with a cut off at 100- 300 keV. The soft 
peak is formed by those photons that escape unscattered 
or have had a small energy change before escape. The 
tail is prominent and extends to high energies even if the 
optical depth is low. It is formed by those few photons 
that picked up a lot of energy from the electrons, mainly 
through head-on collisions. 

Since model fits to the observed X-ray spectra in black- 
hole candidates give blackbody temperatures for the ac- 
cretion disk in the range ~ 0.3 - 1 keV (Tanaka & Lewin 
1995), we have also considered a source of soft photons in 
the disk with blackbody emission of temperature kT = 0.5 
keV. The pointlike source is again placed at the midplane 
of the disk. The emergent spectra are shown in Fig. 3. The 
panels and lines have the same meaning as in Fig. 2. 

The amplitude of the high-energy tail and the cutoff 
photon energy decrease with distance of the photon source 
from the hole. This is expected since the orbital velocities 
encountered by the photons decrease with distance. The 
photon-number spectral index is, however, less dependent 
on the source position. Within statistical errors, it is also 
the same for a monochromatic and a black body input 
spectrum. 

Since the scattering photons sample a velocity range 
that depends on the disk thickness (cf. the discussion in 
section 1.1), the spectra depend not only on the optical 
depth, but also on the actual geometric thickness of the 
disk. This is illustrated in Fig. ^. It shows that the high 
energy tail is strongly suppressed when the aspect ratio 
H/r drops substantially below unity. 

These spectra are reminiscent of the observed spectra 
in black-hole candidates in the high/soft state. The resem- 
blance is strengthened when the photon-number spectral 
index and the ratio of the hard to the soft flux are com- 
pared with observations. 

The photon-number spectral index in the tail varies in 
the range 2-3 and is, within statistical errors, independent 
of the position of the source of soft photons. It depends on 
the vertical optical depth To at r = ri ast ; the thicker the 
medium the softer the spectrum. Figure [| shows the de- 
pendence of the power-law photon-number spectral index 
on the vertical optical depth To . The index was computed 
from the spectra of Fig. |^, in the energy range between 
20 and 100 keV. For low and intermediate optical depths, 
the photon- number index is ~ 2 — 2.5, which is similar to 



those observed in the high/soft state of black-hole candi- 
dates (Tanaka & Shibazaki 1996). 

In addition to the slope of the high-energy tail, we 
can also compare its strength with the observations. We 
do this in terms of a hard-to soft flux ratio. Observations 
of soft state spectra (e.g. Miyamoto et al. 1993; Cui et 
al. 1997), generally give values of 0.01-0.2 for the ratio 
of power law tail to black body fluxes. In Fig. ^ we have 
plotted the 20-200 keV to 1-5 keV flux ratio of our models 
as a function of the vertical optical depth to- 

In a realistic accretion disk, it is expected that the ma- 
jority of the photons in the ultrasoft component come from 
the region of the disk near 6r g . Close to the black-hole 
horizon the vertical optical depth of the disk is relatively 
small and at r 3> 6r g the temperature and emissivity of 
the disk are small. Thus, the squares in Fig. ^ are closer to 
reality than the circles. In other words, the vertical optical 
depth of the disk at the last stable orbit must be greater 
than ~ 3 in order for our proposed mechanism to produce 
spectra similar to the ones observed. 

The energy of the photons escaping to infinity is the 
result of the competition of two processes: the gain of en- 
ergy due to Compton scattering, which increases towards 
the horizon, and the loss of energy due to gravitational 
redshift, which also increases. It can be shown that for a 
thin disk the escaping photon energy cannot exceed m e c 2 . 
Some of the spectra in Fig. || and [| show a cutoff at en- 
ergies a bit above this limit. This can be traced back to 
a simplification made in our assumed velocity field. The 
rotational velocity was taken the same at all z (for the 
same polar coordinate p). On the other hand, in the grav- 
itational redshift the radial coordinate enters spherically, 
i.e. r = \J p 1 + z 2 , which is larger than p. In reality, the 
rotation speed will decrease in step with the redshift. 

In conclusion, we have discussed for the first time the 
Comptonization by orbital motion in an accretion disk 
around a Schwarzschild black hole (as opposed to thermal 
Comptonization in a hot corona or bulk-motion Comp- 
tonization by radially-falling electrons) and have shown 
that high-energy power-law X-ray spectra of the type seen 
in the soft/high states of black hole accreters. 
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